UNITED STATES UTILITY APPLICATION 



HILL & SCHUMACHER 



Title. SUPERRESOLUTION IN PERIODIC DATA STORAGE MEDIA 



Inventors: DARREN KRAEMER Canadian 

33 Himount drive 
Toronto, Ontario 
Canada, M2K 1X3 

BRADLEY SIWICK Canadian 

155 Dalhousie Street 
Apt. 955 
Toronto, Ontario 
Canada M5B 2P7 

R. J. DWAYNE MILLER Canadian 

21 Elmwood Avenue South 
Port Credit, Ontario 
Canada, L5G 3J6 



SUPERRESOLUTION IN PERIODIC DATA STORAGE MEDIA 



CROSS REFERENCE TO RELATED U. S. PATENT APPLICATION 

This patent application relates to United States provisional patent 
application Serial No. 60/258,302 filed on December 28, 2000, entitled 
SUPERRESOLUTION IN PERIODIC OPTICAL DATA STORAGE. 

FIELD OF THE INVENTION 

This invention relates generally to the storage and acquisition of binary 
information that has been stored physically. 

BACKGROUND OF THE INVENTION 

The speed, capacity, and usage of computers have been growing rapidly 
for several decades. With the development of multimedia and digital 
communications, the quantity of data being processed has grown faster than the 
technology that supports it. These factors have given rise to a critical need for 
new modalities of high density memory. For example, data bases at major 
international laboratories such as CERN and medical centers are approaching 
requirements of petabit storage levels per year. Digital media will soon require 
terabit formats, (see S. Esener and M. Kryder, WTEC Panel Report 7 (1999) 
5.2). It may well be that limitations in storing and retrieving data may ultimately 
limit progress in this key sector to science and technology. 

Recently the inventors have investigated the unique properties of a polymer 
photonic crystal with respect to applications as a medium for high-density 3-D 
optical data storage for this expressed purpose, see B. Siwick, O. Kalinina, E 
Kumacheva, R. J. D Miller and J. Noolandi, Appl. Phys. 90 (2001) 5328. This 
new class of materials has the unique property that information can be confined 
in spatially modulated domains with well-defined Fourier components. The 
spatial order imposed on the data storage process is achieved by using 
nanocomposite polymers in which the optically active domains are localized in 
the center of self-assemblying latex particles through a two step (or multi-step) 
growth cycle. Core particles are formed by an optically active dye-labeled 
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monomer. The cores are then subjected to a second stage of polymerization with 
a lower glass temperature that serves as the optically inactive outer buffer. The 
resulting particles show very good uniformity and are essentially monodispersed 
with respect to radii. This property facilitates the formation of hexagonally close 
packed films up to mm thickness with a very high degree of order. Upon 
annealing, the close-packed core-shell particles form a nanostructured material 
with the fluorescent particles periodically embedded into the optically inert matrix 
in a close-packed hexagonal structure, see E. Kumacheva, O. Kalinina, L. Lilge, 
Adv. Mater. 1 1 (1999) 231 . The long range periodic order of these materials is 
demonstrated in Figures 1a), 1b) and 1c). 

To store information in these materials, a two-photon laser scanning 
microscope was used to write information by photobleaching the optically active 
cores and, under much lower fluence, read out the resulting image. The optical 
properties of these nanocomposite polymers are essentially binary in nature, i.e., 
the optically active domains are separated by inactive regions with a square 
modulated cross section determined by the abrupt interface between the high 
and low temperature polymer blends. This binary modulation of the polymer's 
optically properties makes these materials ideal for storing binary information, as 
can be appreciated from Figures 1a), 1b) and 1c). The bleaching of the dye in 
the well defined spatial domains serves as part of the binary code. The sharp 
boundaries and high local concentration of fluorescent dye serves to enhance 
the contrast and readily define the bit storage. Equally important, the writing 
speeds can support GHz address speeds. Overall, a storage density 
approaching 1 Terabit/cc was successfully demonstrated as disclosed in B. 
Siwick, O. Kalinina, E Kumacheva, R. J. D Miller and J. Noolandi, Appl. Phys. 
90 (2001)5328. 

Relative to conventional homogeneous storage media, the nanostructured 
periodic material is shown to increase the effective optical storage density by 
spatially localizing the optically active region and imposing an optically inactive 
barrier to cross-talk between bits. The basic principle that provides this unique 
property is shown schematically in Figure 2. This feature alone leads to a 
significant increase in the storage density of nanocomposite polymers relative to 



homogeneous polymers, (B. Siwick, O. Kalinina, E Kumacheva, R. J. D Miller 
and J. Noolandi, Appl. Phys. 90 (2001) 5328). The aim of the current work is to 
explore the possibility of further taking advantage of the periodic spatial 
properties of these materials. 

Conventional limits in optical resolution (Rayleigh criterion) assume no a 
priori knowledge of the location of an object(s). In the case of periodic materials, 
the objects of interest are confined to lattice points with a well-defined lattice 
spacing such that there is additional information available for signal processing. 
The periodic nature of the signal that would result from such structures is well 
suited to spatial phase sensitive detection methods, i.e. the underlying lattice 
period acts as a reference for signal processing. Signal-to-noise enhancements 
of several orders of magnitude are typically realized in related problems using 
optical heterodyne detection. The question is how far could the spatial resolution 
be improved by explicitly taking advantage of the underlying Fourier spatial 
components to nanocomposite materials. To this end, a post processing 
algorithm was developed that is optimal for extracting binary information from 
images of bit patterns in nanocomposite materials at densities far beyond the 
classic Rayleigh resolution limit. Under realistic noise conditions, this signal 
processing procedure should lead to a density increase of an order of 
magnitude. 

The storage density limit in optical systems is determined by the resolution 
of the imaging system. When observed through an optical instrument such as a 
microscope, the image of an individual bit appears as an intensity distribution. 
The width of this distribution or diffraction pattern is proportional to the 
wavelength X of the light in the imaging system. If the bits are separated by a 
distance X, then there will be little interference between the diffraction patterns of 
the bits. We will refer to this spacing as normal density. As the bit separation is 
reduced, it becomes harder to distinguish individual bits. When the central 
maximum of one bit's image falls on the first minimum of a second bit, the 
images are said to be just resolved. This limiting condition of resolution is known 
as Rayleigh's criterion. In other words, if the distance between bits is smaller 



than 7J2 the microscope will not be able to distinguish between two adjacent bits. 
We will refer to this spacing as the classical limit density. 



(radius 4cm) has a density of 13 Mbytes/cm 2 . For such systems which use a 780 
nm (infrared) laser the normal density is around 20 Mbytes/cm 2 . Hence, 
commercial CD-ROMs operate at about 60% of the normal density. The bit 
density of a DVD that uses a 650nm (red) laser is 7.8x1 0 8 bits/cm 2 . For this 
wavelength the normal density is around 1 .64x10® bits/cm 2 , and the classical limit 
density, is 9.5x1 0 8 bits/cm 2 . The commercial DVD's operate near 80% of the 
classical limit density. Although the DVD's density is a great improvement (~7 
times more data per disc) over the older CD-ROM technology, it is still limited by 
diffraction as expressed in the Rayleigh criterion.. 

When two diffraction patterns are superimposed, in order to satisfy the 
Rayleigh criterion, they must be separated by more than half the width of the 
central peak. 



where z is x, R is y, and Xq is z. The minimum resolved separation distance 
(nj) is therefore proportional to half the wavelength. This relation defines the 
maximum density that can be stored in the image plane and is responsible for 
the r 2 scaling for 2-d optical storage. By taking advantage of confocal 



above relation then defines the storage density for one page in the optical 
readout for 3-d memory and the confocal parameter defines the density of 
these pages along the optic axis. The density along the optic axis can be 
increased to a certain extent by using theta scans (see Steffen Lindek, Rainer 
Pick, Ernst H K Stelzer Rev. Sci. Instrum. 65, 3367 (1994)) however, the 
storage density scales quadratically for in-plane components and provides far 
greater gain in storage density. 

It would be very advantageous to have a technique that allows data to 
be read at resolutions beyond the Rayleigh criterion. 



A CD-ROM that holds a maximum of 650 Mbytes over an area of 50cm 2 




[1] 



imaging to provide a 3 dimension, the maximum density scales as X . The 



SUMMARY OF THE INVENTION 

Accordingly, one object of the present invention is to provide a 
technique which allows data to be read from a periodic data storage media at 
resolutions beyond the Rayleigh criterion. 

Another object of the present invention is to provide a method to 
resolve structures in images beyond the Rayleigh criterion, as well as to 
extract the binary information contained therein. Error Free determination of 
this information at bit densities exceeding the classical limit would thus imply, 
an effective "super-resolution". 

For a periodic lattice memory system, the superresolution or MMD 
method disclosed herein is achieved by post-processing the images of two- 
dimensional pages of memory, either recorded simultaneously or time- 
averaged (as with one-dimensional beam systems). This post-processing 
technique, described herein, is referred to as matrix-method deconvolution, 
since it compensates for the effects of the imaging system's point spread 
function. Prior knowledge of a system's point spread function and interbit 
spacing is used. Also the center co-ordinate of at least one bit must be 
determined to set the effective phase of the spatially modulated signal. 

The present invention provides a method of reading binary information 
stored in a storage medium, comprising: 

a) providing a storage medium having n memory-centers each with a 
known position and the memory-centers having substantially the same 
physical dimensions; 

b) accessing said storage medium with an addressing system and 
measuring for each memory-center a scalar signal intensity l m emitted from a 
pre-selected region which is centered on the known position of said memory- 
center; and 

c) extracting the stored binary information by calculating bit values b n 
for all memory-centers using an equation B = C" 1 1 /l 0 , wherein l 0 is a 
predetermined normalizing factor, i = (h, l 2 ,...,l n ) is an array of said scalar 
intensities f6r all memory-centers, and B = (bi, b 2 , b n ) is an array of bit 



values, and C is a predetermined cross-talk matrix of n 2 elements where each 
element represents a cross-talk between said pre-selected regions. 

In this aspect of the method, the value of each matrix element may be 
defined as a function of a spacing between memory-centers i and j given by 

Cij= f(f) = f(| £i - = f(Rij) 
where f(r') is defined as a cross-talk function, and wherein said cross-talk 
matrix C is calculated by applying said cross-talk function to each element of a 
matrix R that contains all inter-memory-center spacings R i j=r' = |rs - rj|. 

In this aspect of the method the cross-talk function f(r') may be derived 
from an intensity distribution within a preselected region l 0 (flm), 

f{f) = ]dql,{q)^l,{q) 

where q m defines coordinates of the intensity distribution of the preselected 
region of the mth memory-center. 

In this aspect of the invention a binary value for each memory-center is 
calculated from a corresponding bit value by a process wherein the ni highest 
bit values are assigned a binary value of T and all others are assigned a 
binary value '0' based upon an equation relating the population of '1' valued 
memory-centers, 

N r. jtotal 

Plo I 0 

In another aspect of the invention there is provided a method of reading 
binary information stored in a storage medium, comprising: 

a) providing a storage medium having n memory-centers each with a 
known position and the memory-centers having substantially the same 
physical dimensions; 

b) accessing said storage medium with an addressing system and 
measuring for each memory-center a scalar signal intensity l m emitted from a 
pre-selected region which is centered on the known position of said memory- 
center and having an intensity distribution defined by an impulse response of 
the addressing system and an effective distribution of the signal stored within 



the addressed memory-center; and 

c) extracting the stored binary information by calculating bit values b n for 
all memory-centers using an equation B = C 1 i /l 0 , wherein l 0 is a 
predetermined normalizing factor, i = (h, l2,--.Jn) is an array of said scalar 
intensities for all memory-centers, and B = (bi, b 2 , b n ) is an array of bit 
values, and C is a predetermined cross-talk matrix of n 2 elements where each 
element represents a cross-talk between said pre-selected regions. 

The present invention also provides a method of reading binary 
information stored in an optical storage medium, comprising: 

a) providing an optical storage medium having n memory-centers each 
with a known position and the memory-centers having substantially the same 
physical dimensions; 

b) accessing said optical storage medium with an optical addressing 
system and measuring for each memory-center a total optical intensity l m 
emitted from a pre-selected region which is centered on the known position of 
said memory-center and having an optical intensity distribution within a single 
pre-selected region l 0 (q) defined by a point spread function of the optical 
addressing system and an intensity distribution of the memory-center itself 
defined by an optical response of a single memory-center as imaged through 
an idealized optical addressing system having an infinitely small point spread 
function; and 

c) extracting the stored binary information by calculating bit values b n for 
all memory-centers using an equation B = C" 1 1 /l 0 , wherein l 0 is a 

predetermined normalizing factor, I = (h, l2,.-.Jm In) is an array of said scalar 

intensities for all memory-centers, and B = (bi, b 2 , b n ) is an array of bit 
values, and C is a predetermined cross-talk matrix of n 2 elements where each 
element represents a cross-talk between said pre-selected regions. 

An advantage of the method of the present invention is that it removes 
the data storage limit classically imposed by the Rayleigh criterion with the 
only limit of data retrieval being related to noise conditions rather than 
density. 



BRIEF DESCRIPTION OF THE DRAWINGS 

The method of the present invention will now be described, reference 
being had to the accompanying drawings, in which: 

Figure 1a shows a structure of the addressed plane of a 
nanocomposite memory storage medium located at the distance 200 microns 
from the top surface, the scale bar is 5 microns and the particles have 500 
nm diameter cores and 250 nm thick shells; 

Figure 1b) shows the structure of the addressed and adjacent layers of 
the nanocomposite storage medium of Figure 1a) studied by 2-photon 
confocal microscopy; 

Figure 1c) shows fluorescence intensity profiles for the addressed and 
neighboring planes of Figures 1a) and 1b); 

Figure 2 shows how a memory medium made of nanostructured 
materials significantly reduces the cross-talk in the writing and reading 
processes by spatial isolation/separation of the active memory centers 
compared to homogenous memory storage materials; 

Figure 3a shows the Regions of Interest (ROI) with no overlap 
between ROI's 

Figure 3b(i) shows the memory-centers or bits are close enough that 
ROI's overlap a little, 3b(ii) shows the resulting cross-talk between the ROI 
from the central ON bit; 

Figure 3c(i) shows the memory centers close enough so there is high 
overlap of the ROI's, 3c(ii) shows the resulting cross-talk between the ROIs in 
3c(i); 

Figure 4a illustrates the definitions of r, the spacing between the 
memory centers or bits, and R, the radius of the diffraction pattern; 

Figure 4b shows the basis vectors and the lattice translations for the 
periodic storage memory where the dark circles represent the array of 
periodic memory centers; 

Figure 5 (top left ) shows a laser confocal fluorescent microscopy 
image of a nanocomposite particle array in which the memory centers have a 
core diameter 650 ± 20 nm and shell thickness 200 ± 5nm . A, F i U orescnce~500 
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nm, Xtwo- P hoton= 844nm, the resolution is approximately 256x256 samples, in 
which a data pattern has been photo-bleached into material (top right); after 
filtering and deconvolution approximate Gaussian point spread function is 
shown in the bottom right, and at the bottom left a simulation of equivalent 
5 data with a sine squared basis bit is shown with simulated point-spread 
function of diameter, 750 nm and signal to noise ratio of 10; 

Figures 6a, 6b, 6c and 6d show increasing overlap in random data at 
various densities; 

Figure 7 shows bit distributions, when the MMD method is applied to 
10 simulated images the output values of the B n vector tended to take on two 
distributions around 1 and 0; and 
O Figure 8 shows plots of cross-talk % versus the signal to noise ratio for 

2 either simulation show the regions in which the presen MMD method of the 

present invention can improve the resolving power of memory acquisition 

y 15 systems. 

W 

£ DETAILED DESCRIPTION OF THE INVENTION 

H A goal of the present method is not to resolve structures in the images 

n beyond the Rayleigh criterion, but rather to extract the binary information 

H s 20 contained therein. Error Free determination of this information at bit densities 
exceeding the classical limit would thus imply, an effective "super-resolution" . 

In the present method, for a periodic lattice memory system, 
superresolution is achieved by post-processing the images of two- 
dimensional pages of memory, either recorded simultaneously or time- 
25 averaged (as with one-dimensional beam systems). This post-processing 
technique, described herein is referred to as matrix-method deconvolution 
(MMD), since it compensates for the effects of the imaging system's point 
spread function. Prior knowledge of a system's point spread function and 
inter-bit spacing is needed. Also the center co-ordinate of at least one bit 
30 must be determined to set the effective phase of the spatially modulated 
signal. 
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Matrix-Method Deconvolution (MMD) 

We define the Region of Interest (ROI or pre-determined region) in an 
image of a bit to be that area encompassing the first (and most significant) 
maximum of the bit pattern. Figure 3a shows the situation of no overlap 
between ROI's, Figure 3b(i) shows memory centers being close enough so 
that ROI's overlap a little and Figure 3b(ii) shows the resulting cross-talk 
between the ROI's from the central ON bit or memory center. Figure 3c(l) 
shows highly overlapping ROI's and Figure 3c(ii) shows the resulting cross- 
talk between the ROI's. Since this information can be reconstructed as a 
Fourier series of the lattice period, the pattern appears as an effective 
diffraction pattern to this periodic function and will be referred to as the bit 
diffraction pattern. The total detected intensity l m within the mth ROI ultimately 
determines the corresponding binary value. At densities in which there is an 
overlap (cross-talk) between the ROI's, the total intensities l m will increase by 
an amount related to the proximity between bits, such that: 

n 

where B N is the binary value of the nth bit, l 0 is the maximum intensity within 
an isolated ROI and c xn = c nx = / (r n -r x ) is the cross-talk coefficient. 

If we include all of the bits in an array of N bits, we have a system of N 
equations and N unknowns. 

h = l 0 x (Bi + c 12 B 2 + C13B3 ...+Ci N B N ) 
l 2 = lo x (C21B1 + B 2 +C23B3 ...+c 2 nBn) 

1 
i 
I 
I 

In = lo x (CN1B1 + CN2B2 + c N3 B 3 ...+B N ) [3] 
In matrix form, these equations can be expressed as 

C nm B m = l n / lo [4] 

where C is an n x n matrix that shall be called the cross-talk matrix. The 
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overlap leading to cross-talk between bits is shown schematically in Figure 3. 

To extract data from an image, we measure the intensities within each 
ROI, and calculate the bit values using the inverse of the previous equation: 

B m =C nm 1 L/lo [5] 
C and l 0 are constant for a given combination of readout frequency and bit 
density, and can be calculated a priori if the function / (r jk ) is applied to an 
array of inter-bit distances, n k =n - r k . 

Matrix Construction and the Cross-talk Function 

The cross-talk function /(r !k ) is determined by the intensity distribution 
l 0 (x,y) of a bits image within its ROI, see Figure 4a, 

= 4 {/2 f ( x)l o( x >y) d y dx [6] 

or for discrete pixels 

/(>•) = 4X§/.(x,jO 17] 

r/2 0 

With higher resolutions, it was possible to approximate the cross-talk function 
/(ry) with a polynomial fit to speed up calculations, but for lower resolutions 
the pixels make this smooth function very jagged, and instead, a look-up 
table was generated. 

For each unique combination of bit size, spacing and point-spread 
function, a different matrix has to be generated. Since there are N 2 
components in Cy, the theoretical matrix construction is a second order 
algorithm. However, it can be reduced to an Nlog N order algorithm by taking 
advantage of the periodic lattice. To generate the cross-talk matrix, an array 
was created for each bit containing the distance between it and the other bits; 

=|rjj| is also an nxn matrix. Then the cross-talk function /(nj) was applied to 
each element of R. 

The first row of R corresponds to the distances between the first bit 
and all others bit. This is calculated using integer combinations of lattice 
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basis vector as shown schematically in Figure 4b, 

Rij = A s a + Bjb + CjC. [8] 
Then the jth row was calculated by transforming the origin from the first bit to 
the jth, i.e., 

|R jk | = |Rik-Rijl [9] 
Once density specifications are chosen for a system, the cross-talk matrix 
can be experimentally determined by recording the ROI intensities when only 
one bit in the system is "on.". The intensities measured within all other ROI 
correspond to one row of the cross-over matrix. All subsequent rows are 
lattice vector translations of the first. 

As stated above, the cross-talk matrix has N 2 elements. If N 2 
instructions must be executed to access N bits, the speed of the memory 
system will be very slow. However, for values of N much larger than the 
number of overlapping diffraction patterns within a single ROI, the overlap 
matrix will be very sparse. As the size of N increases the matrix becomes 
increasingly sparse and fewer calculations are necessary to extract the 
information. 

It should be noted that when memory is within the classical limit 
density, it is possible to extract the bits with a simple algorithm that compares 
the intensities of the ROI to a set threshold of say nl 0 /2 where n is the number 
of nearest neighbors. Yet even at these classical densities, there are 
advantages to using the MMD extraction procedure. It calculates each bit as 
a function of the intensities of the surrounding bits, like a moving 
average/autogain algorithm, and further enhances the signal to noise as will 
be discussed below. 

Simulation of Experimental Results 

To simulate the two-photon confocal flourescence microscopy images of 
data stored in nano-structured polymer, we assumed a bit basis function: the 
first peak of a radially symmetric cosine squared function, 
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S(x, y) = Cos 1 ^- 4(x -«A) 2 +(y- %f ) for (x-«/ 2 ) 2 + (y-«A) 2 <R 2 
2R 

= 0 (ignores second order maxima) for (x - R / 2 ) 2 + (y - R A) 2 > R 2 

[10] 

This function was used as the basis to construct simulated close packed 
hexagonal lattices of random binary data and is shown in Figure 5. 
Simulations were conducted by two different methods, to demonstrate 
different situations. In both methods, the diffraction of the bit basis was 
calculated by convolving with an approximate Gaussian point-spread 
function. 

I 0 (x 9 y) = I PSF (x 9 y)*S (x 9 y)= jjl psp (x 9 y)S {x'-x 9 y'-y)dx'dy' 

where I PS f is a spatial distribution of the impulse response and S(x,y) is the 
optical intensity distribution of the single memory-center as measured through 
an idealized optical addressing system. 

Simulation A 

The properties of the imaging system were held constant and the 
information storage density was varied. Hence, the radius of the diffraction 
pattern was kept constant and the bits were moved closer together. The 
samples per unit area, the area of the ROI, and the magnitude of each 
diffraction pattern is held constant. When the profile of diffraction patterns are 
overlapping by 50% , the distance between nearest neighbor bits is equal to 
half the diameter of the effective diffraction pattern (A/2) , and it becomes 
difficult to visually distinguish nearest neighbor bits. This point is illustrated 
schematically in the series showing increasing overlap in Figures 6a, 6b, 6c 
and 6d. This is the classical limit density. At 75% , the distance between bits 
is a quarter the diameter of the diffraction pattern (AY4). An image of the same 
quantity of data is now one sixteenth the size of the uncompressed lattice. 
This translates to a 3-d density improvement over normal density of 256 
times. 
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Simulation B 

The density of stored information was held constant arid the properties 
of the imaging system (PSF) were varied. Alternatively this procedure could 
be applied to a given imaging system, a given data density, by varying the 
wavelengths of light Xq. The distance between the bits is held constant; while 
the width of the diffraction pattern is changed. As the point-spread function 
widens, the maximum intensity goes down. When the radius of the diffraction 
pattern is equal to the spacing between the bits, the effective cross-talk is 
50%. This is the classical limit density for such a system,A=2A,o. At 75% 
effective cross-talk, X=3X 0 , the diameter of the diffraction pattern is 3 times 
the distance between bits or, 

d 

Effective Cross-talk = (1- X 0 fX) *100 = (1 ) *100 [11] 

2 G PSF 

where d CO re is the core diameter and G PS f is the FWHM of a gaussian point 
spread function. It should be emphasized that for this situation, an effective 
cross-talk of 100% requires an infinite point spread function; whereas in the 
previous simulation, 100% cross-talk implies that both bits were in the same 
location. 

When the MMD was applied to simulated images the output values of 
the B n vector tended to take on two distributions around 1 and 0, as shown in 
the bit distributions of Figure 7 which is a series of MMD output histograms. 
As the noise was increased, the widthof both ON and OFF peaks in the 
histograms become wider, until a noise threshold was reached, at which point 
the two bands of values begin to overalp and the algorithm is no longer 100% 
accurate. Higher resolutions are affected less by the noise because it is 
averaged whe nthe intensities Ij are calculated from ROI with more bits 
contained within them. To extract information from the output, binary values 
must be assigned to each bit. Note that the population of 1's and 0's, is 
determined from the total intensity within the image. 

n,=tlTlL <h-N-n, ' [12] 



14 



The ni highest valued bits are T and the remaining n 2 =N-ni bits are 
considered '0'. 

The resolving power of the MMD method was tested in both situations 
by adding different kinds and levels of noise to the simulated images. The 
signal-to-noise ratio was defined as the ratio between the signal power (total 
intensity) and the noise power. The results summarized in Table 1 were 
obtained by simulating arrays of 256 bits with each bit having a diameter of 
16 pixels. These values were held constant throughout the simulations, so as 
not to skew the results. The effects of changing the numbers of pixels per bit 
disappear for large enough numbers of bits. The speed of the simulations 
was increased by choosing sampling parameters that were powers of two 
(allows an FFT algorithm to be used). 

The noise threshold for error free determination of various 
configurations of both simulation A and B were measured by slowly 
decreasing the amount of noise until no errors were detected. The regions of 
Error Free Signal (EFS) demonstrated with the MMD are shown in Figure 8. 

Table. 1 Summary of Noise Threshold Results: 
Simulation A 



Threshold SNR [dB] 



2D 
Density 
Improvem 
ent over 
normal 
density 



Effective 
Cross-talk 



White 
Noise 



[%] 



p/ P normal 



o 

30 
50 
70 
75 
90 
99.9 



-3 
-3 
-3 
20 
30.5 
46. 

[00] 
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2 

4 

11 

16 
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Simulation B 



Effective 

Cross-talk 

[%] 



Threshold SNR [dB] 
White 
Noise 



2D 
Density 
Improvem 
ent over 
normal 
density 



0 

30 
50 
67 
75 
80 
83 
99.9 



-3 
-3 
-3 
-3 
10 
27.5 
49 

[CO] 



0 

2 

4 

9 

16 

25 

36 

M 



The calculated improvements in information density are rather 
remarkable. This potential increase in bit density that can be accessed 
through such a binary extraction method should be gauged relative to other 
approaches being explored. A great deal of technical effort is presently being 
expended towards the development of blue diode lasers to increase optical 
data storage by approximately a factor of two. By simply changing the 
material to a nanocomposite and using the MMD approach it is possible to 
increase the density substantially beyond this target. The effectiveness of this 
technique is due to the following factors; knowledge of the effects of the 
imaging system's point spread functions, and the structure of the storage 
materials, i.e. the spatial location of each bit in the images. 

It will be appreciated that the increase in data storage density is 
determined by the signal-to-noise in the read out. In particular, the width of 
the distributions in the algorithms output is related to the amount of noise in 
the image. This shows that the maximum resolvability of MMD will be 
dependant upon the system's signal-to-noise ratio, not the diffraction limit 
defined solely by the optics. The noise effects may be reduced with spatial 
and temporal filters to further suppress background noise and to build in 
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redundancy into an error correction process. 

There are several significant advantages of the matrix-method 
deconvolution (MMD) method disclosed herein. It can be exploited in a 
number of different applications. It can be used to read out data at longer 
wavelengths than the writing process. In order to write information beyond 
the diffraction limit of the reading wavelength requires that the data be written 
at shorter wavelengths as this step is still defined by diffraction limits of the 
optical imaging system. The writing process can be done with more 
expensive "blue lasers" to encode higher density information and MMD 
enables the read out at longer wavelengths with more cost effective systems 
amenable to multiple users. The MMD method can be used for reading low 
density memory, in a system that has high noise levels. It can be used to 
read from storage media systems with limited sampling capability. The MMD 
is only dependant on the total intensity within a ROI, hence the minimum 
sampling required is one sample per ROI. The matrix method of 
superresolution can work in either one, two, or three dimensions, whether the 
pre-selected regions of interest are lines, circles or spheres. 

To summarize, a method has been disclosed hereinto achieve error- 
free retrieval of binary information an order of magnitude beyond the Rayleigh 
limit. The method disclosed herein to achieve substantially error-free retrieval 
of binary information an order of magnitude beyond the Rayleigh limit, known 
as the matrix-method doconvolution method (MMD), was tested on images 
that simulated an existing optical memory system. It was found that the MMD 
method successfully compensatesfor the limitations of diffraction. Within 
current signal to noise limits in data read out, increases in optical density of 
more than 64 over the classical density are achievable, without considering 
further signal conditioning and optimization. 

The exploitation of the intrinsically periodic nature of nanocomposite 
materials gives memory storage media made of these nanocomposites an 
advantage over homogeneous memory storage media in this regard. 
Notwithstanding, it will be understood that the present method can be used 
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with homogenous or non nanocomposite media as well, although the biggest 
improvement is expected using nanocomposite memory storage media. 

The underlying period of the nanocomposite lattice provides an 
effective spatial reference for signal processing and extracting additional 
information. The combination of novel periodically photoactive materials and 
image processing may provide the solution to our ever-increasing demand for 
more memory at higher densities. Further, nanocomposite materials bring an 
extra dimension for optimizing the photoactive processes for a variety of 
applications. By spatially varying the viscosity and polarity of the 
nanocomposite, it is possible to increase photochemical quantum yields, 
create enhanced photo-electric index modulation, and minimize volume 
contraction in the recording process. All of these attributes can be further 
enhanced with respect to information density by explicitly exploiting the 
periodic nature of these materials, as illustrated with MMD approach for the 
simplest example of binary code. 

One type of nanocomposite material is the periodic array of nano- 
particles being a polymer matrix comprising a three dimensional array of rigid 
polymeric cores embedded in a substantially transparent shell-forming 
polymer, each memory-center comprises a photosensitive constituent 
associated with each nano-particle. The photosensitive constituent may be 
chromophores on, within or otherwise associated with the cores. More 
particularly the chromophores may be fluorescent molecules. The information 
is stored optically within the memory-centers in the storage medium and the 
addressing is an optical addressing system. 

While the method has been described in detail with respect to optical 
storage media and addressing systems, it will be appreciated that the method 
of the present invention is not restricted to optical storage and addressing 
systems. For example, the storage medium and addressing system may be a 
magnetically based system in which the information is stored magnetically 
within the memory-centers in the storage medium, and wherein l m is the 
magnetic intensity within the pre-selected region of the mth memory-center. 
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In this case each memory-center comprises a magneto-sensitive constituent 
associated with each nano-particle when the storage medium is comprised of 
a periodic array of nanoparticles. The polymeric cores may be rigid polymeric 
cores such as, but not limited to, latex spheres. 

Thus in the broadest aspect the present invention provides a method 
of reading binary information stored in a storage medium, comprising 

providing a storage medium having n memory-centers each with a 
known position and the memory-centers having substantially the same 
physical dimensions; 

accessing the storage medium with an addressing system and 
measuring for each memory-center a scalar signal intensity l m emitted from a 
pre-selected region which is centered on the known position of the memory- 
center; and 

extracting .the stored binary information by calculating bit values b n for 
all memory-centers using an equation B = C" 1 I /l 0 , wherein l 0 is a 
predetermined normalizing factor, i = (h, 12,-Jn) is an array of the scalar 
intensities for all memory-centers, and B = (bi, b 2) b n ) is an array of bit 
values, and C is a predetermined cross-talk matrix of n 2 elements where each 
element represents a cross-talk between the pre-selected regions. 

The value of each matrix element may be defined as a function of a 
spacing between memory-centers i and j given by 

C i rf(0 = f(|r i -r j |) = f(Rij) 
where f(r') is defined as a cross-talk function, and wherein the cross-talk matrix 
C is calculated by applying the cross-talk function to each element of a matrix 
R that contains all inter-memory-center spacings Rjj=r' = |ij - rj|. 

The cross-talk function f(r') may be derived from an intensity 
distribution within a preselected region lo^m), 

/(rO = ^/ 0 (^)n/ 0 %) 

where q m defines coordinates of the intensity distribution of the preselected 
region of the mth memory-center. 

A binary value for each memory-center is calculated from a 
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corresponding bit value by a process wherein the ni highest bit values are 
assigned a binary value of '1' and all others are assigned a binary value '0' 
based upon an equation relating the population of '1' valued memory-centers, 

N TJ rtotal 

jiio i o 

In all embodiments of the storage medium and addressing system, the 
storage medium may be a 1-, 2- or 3-dimensional storage medium. The 
storage medium may be addressed in 1-, 2- or 3-dimensions using the 
addressing system. 

The storage medium may not necessarily be a nanocomposite based 
storage medium and instead my be a homogeneous optical storage medium in 
which the memory-centers are "written-in" with pre-selected dimensions. 

Similarly, the storage medium may be a homogeneous magnetic 
storage material. 

The foregoing description of the preferred embodiments of the invention 
has been presented to illustrate the principles of the invention and not to limit 
the invention to the particular embodiment illustrated. It is intended that the 
scope of the invention be defined by all of the embodiments encompassed 
within the following claims and their equivalents. 
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